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ABSTRACT 

The Angstrom Project is using a distributed network of two-meter class telescopes to conduct a high cadence 
pixel-lensing survey of the bulge of the Andromeda Galaxy (M31). With the expansion of global telescope 
network, the detection efficiency of pixel-lensing surveys is rapidly improving. In this paper, we estimate the 
detection rate of binary lens events expected from high-cadence pixel-lensing surveys toward M31 such as 
the Angstrom Project based on detailed simulation of events and application of realistic observational condi- 
tions. Under the conservative detection criteria that only high signal-to-noise caustic-crossing events with long 
enough durations between caustic crossings can be firmly identified as binary lens events, we estimate that the 
rate would be T\, ~ (7 — 15)/t,(iV/50) per season, where fi, is the fraction of binaries with projected separations 
of 10~ 3 AU < d < 10 3 AU out of all lenses and N is the rate of stellar pixel-lensing events. We find that detected 
binaries would have mass ratios distributed over a wide range of q > 0.1 but with separations populated within 
a narrow range of 1 AU < d < 5 AU. Implementation of an alert system and subsequent follow-up observations 
would be important not only for the increase of the binary lens event rate but also for the characterization of 
lens matter. 

Subject headings: gravitational lensing - galaxies: individual (M31) 



1, INTRODUCTION 

Surveys to detect transient variations of stellar brightness 
caused by gravitation microlensing have been and are being 
conducted towards various star fields. These fields in clude the 
Magellanic Clouds (MACHO: lAlcock et all d2000h. EROS: 
Afonso et al. (2003)) , Galactic bu l ge (M ACHO: lAlcock et alJ 
d2001|) EROsTfHamadache et al. (2006), OGLE: ISumi et all 
( 2006). MO A: iBond et al. n200lb), a nd M31 (AGAPE: 
lAnsari et all (119991). POI NT- AGAPE: ICalchi Novati et all 



d2005h. VATT-Colo mbia: lUglesich et al.1 (120041). M EGA: 
Ide Jong et alJ d2004h . WeCAPP: iRiffeser et al.l (120031) ). The 
number of lensing events detected so far from these surveys 
is about 3000. Most of these events were detected toward 
the Galactic bulge field and the number of events detected 
toward the Magellanic Cloud and M31 fields, which are sev- 
eral dozens toward the individual fields, have meager contri- 
butions to the total number of events. Towards the Magellanic 
Cloud field, the line of sight passes mainly through the halo 
of our Galaxy and thus the low detection rate is attributed to 
the scarci ty of massive compact halo objects that can work 
as lenses (lAlcock et al.ll2000t iTisserand. etaT]|2006l) . On the 

1 Department of Physics, Institute for Basic Science Research, Chungbuk 
National University, Chongju 361-763, Korea 

2 Astrophysics Research Institute, Liverpool John Moores University, 
Twelve Quays House, Birkenhead, Merseyside CH41 1LD, United Kingdom 

3 Jordell Bank Observatory, University of Manchester, Macclesfield, 
Cheshire SKI 1 9DL, United Kingdom 

4 Korea Astronomy and Space Science Institute, Hwaam-Dong, Yuseong- 
Gu, Daejeon 305-348, Korea 

5 Department of Astronomy, The Ohio State University, 140 West 18th 
Avenue, Columbus, OH 43210 

6 Department of Physics and Astronomy, Seoul National University, Seoul 
151-742, Korea 

7 Ulugh Beg Astronomical Institute, Uzbek Academy of Sciences, 
Tashkent, Uzbekistan 

8 Corresponding author 



other hand, the line of sight toward M3 1 passes through the 
dense stellar region of M3 1 and thus the small number of de- 
tected stellar lensing events is mainly due to the low detection 
efficiency of the surveys. 

With the expansion of global telescope network, however, 
the detection efficiency of M3 1 lensing surveys is expected 
to greatly improve. For example, a new M31 pixel-lensing 
survey, the Andromeda Galaxy Stellar Robotic Microlensing 
(Angstrom) Project 9 will be able to achieve a monitoring fre- 
quency of ~ 5 observations per 24-hour period by using a net- 
work of telescopes, including the robotic 2m Liverpool Tele- 
scope at La Palma, Faulkes Telescope North in Hawaii, 1.8m 
telescope at the Bohyunsan Observatory in Korea, the 2.4m 
Hiltner Telescope at the MDM Observatory in Arizona, and 
1 .5m at the Maidanak Observatory in Uzbekistan. Intensive 
monitoring programs such as Angstrom a re expected to detec t 
events with a rate of ~ 100 per season (iKerins et al.1 12006). 
Among them, a considerable fr action would be caused by bi- 
naries. iBaltz & Gond olo (2001) pointed out that the rate of 
binary events relative to single lens events is higher in pixel- 
lensing because many of binary events involve caustic cross- 
ings during which the magnification is high and thus more 
detectable. In this work, we estimate the detection rate of bi- 
nary lens events expected from a high cadence pixel-lensing 
survey toward M3 1 based on detailed simulation of events and 
application of realistic observational conditions. 

The paper is organized as follows. In § 2, we describe ba- 
sics of binary lensing. In § 3, we investigate the dependency 
of the detection efficiency on the binary lens parameters. In § 
4, we estimate the event rate based on the simulation of M3 1 
pixel-lensing events. We describe the details of the simulation 
and the criteria applied for the selection of binary events. We 
also investigate the characteristics of the binaries detectable 

9 http://www.astro.livjm.ac.uk/angstrom/ 
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from the survey. We end with a brief summary of the results 
and discussion in § 5. 

2. BASICS OF BINARY LENSING 

General relativity predicts that a light ray passing by a stel- 
lar object is deflected. If a source star is gravitationally lensed 
by a binary lens, the equation of lens mapping from the lens 
plane to the source plane is expressed as 



7=1 



rrij/M 

Z-Zl/ 



(1) 



where £ = £ + if], zlj = xlj + tyLj, and z=x+iy are the com- 
plex notations of the source, lens, and image positions, re- 
spectively, z denotes the complex conjugate of z, nij are the 
masses of the individual lens componen ts, and M = m\ +rri2 is 
the total mass of the system dWitti ll990). Here all lengths are 
normalized to the radius of the Einstein ring of the total mass 
of the system. The Einstein radius is related to the physical 
parameters of the lens system by 



fAGM\ ilA f\ 1_ 



1/2 



(2) 



where and Ds are the distances to the lens and source, 
respectively. Due to lensing, the source star image is split 
into several segments. For binary lensing, the number of im- 
ages, Nj, is either three or five depending on the source po- 
sition with respect to the lens position. The lensing process 
conserves the source surface brightness and thus the magni- 
fications A, of the individual images correspond to the ratios 
between the areas of the images and source. For an infinitesi- 
mally small source element, the magnification is 



A,= 



dz dz 



(3) 



Then, the total magnification corresponds to the sum over all 

images, A = Ylf'^i- 

The main new feature of binary lensing compared to sin- 
gle lensing is the formation of caustics. Caustics are the set 
of positions in the source plane on which the magnification 
of a point-source event is infinite. The set of caustics form 
closed curves, which are composed of multiple concave line 
segments (fold) that meet at points (cusp). When the source 
enters the caustic curve, two new images appear and the num- 
ber of images changes from three into five. Due to the diver- 
gent nature of the magnification near a caustic, the light curve 
during the caustic crossing of a binary lens event is character- 
ized by a sharp spike. Since the caustic curve is closed, the 
number of caustic crossings is a multiple of two. The light 
curve between a set of two caustic crossings is characterized 
by its distinctive "U"-shape. 

The number and shape of caustics vary depending on the 
separation, s (normalized by #e), and mass ratio, q, between 
the two lens components. If the separation is substantially 
smaller than the Einstein radius, s -C 1, there exist three sets 
of caustics. One big caustic is located close to the center of 
mass of the binary and the other two tiny ones are located 
away from the center of mass on the heavier lens side. If 
the separation is substantially larger than the Einstein radius, 
i> 1, on the other hand, there exist two sets of caustics. They 
are located close to the positions of the individual lens com- 
ponents but slightly shifted toward the lens component of the 



nij/M 
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other side. The amount of the shift is 

Az L ., = sgn(z w -z L ,,) 

where zlj is the position of the lens component on the same 
side of each set of caustics while mj/M and zlj represent 
the mass fraction and position of the oth er lens component 
locate d on the opposite side, respectively dDi Stefano & Maol 
1996). If the binary separation is equivalent to the Einstein 
ring, s ~ 1.0, the caustic curve forms a single large closed 
figure with its center roughly at the center of mass of the bi- 
nary. In FigureQ] we present the caustic pattern for various bi- 
nary lenses with different separations. The dependence of the 
caustic size on the mass ratio is weak. As a result, the caustic 
size is not negligible even for an extreme case of a star-planet 
binary lens pair with q < 10~ 3 . Thi s makes microlensing a n 
efficient method for planet searches (Mao & Pacz vriskill991h . 

3. DETECTION EFFICIENCY 

Due to the variation of the caustic characteristics depend- 
ing on the binary lens parameters, the detection efficiency of 
binary lens events also depends on these parameters. In the 
section, we investigate the dependencies of the detection effi- 
ciency on the binary lens parameters. 

We estimate the efficiency for an example M3 1 event. We 
choose a representative event based on the simulation of M3 1 
pixel-lensing events conducted by lKerins et alJ ((2006). In the 
simulation, M3 1 matter density distribution is based on stan- 
dard double-exponential disk plus triaxial bulge models. The 
velocity is modeled by an isotropic Maxwellian distribution. 
The mass function of lenses are modeled by a broken power- 
law represented by 

' n(m/0.5 MoT 14 , for 0.08 M <m<O.5M , 
K(ffi/O.5M )" 235 , for 0.5 M Q <m< 1.0 M , 



<1>(ot) : 



(5) 

where n is a proportional constant. According to this simula- 
tion, the total number of events expected to be detected by the 
Angstrom survey is ~ 30-60 events per season over an area 
of 1 1 x 1 1 arcmin 2 of the M3 1 bulge field under the assump- 
tion that all lenses are composed of a single stellar compo- 
nent. Based on these events, we select a representative event 
as the one caused by a binary lens system with a primary lens 
mass of m = 0.4 Mq. The assumed distances to the source 
star and lens are Ds = 780 kpc and = (780-0.8) kpc, re- 
spectively. For a representative value of the Einstein timescale 
(time needed for the source to transit the Einstein radius of the 
pr imary lens), we ado pt = 10 days based on the simulation 
of iKerins et alJ (;2006). For the source star, we choose a giant 
with a brightness / ~ 24.4. We consider finite-source effects 
by adopting a source radius of = 10.0 Rq. We assume an 
/-band background brightness of /i/ = 17.6 mag/arcsec 2 . We 
let the mass ratio and separation between binary lens compo- 
nents vary as free parameters. 

Once the lensing parameters are set, we then produce lens- 
ing events. Due to the high background flux level, detec- 
tions of lensing events toward M3 1 field are limited only to 
those with high magnifications. In the background-dominated 
regime such as the M31 field, the signal-to-noise ratio of a 
pixel-lensing event is 



S/N-- 



F s tUp(A-l) 

cV2 ' 



(6) 



where Fs and Fa represent the flux from the source and back- 
ground, respectively, and f exp represents the exposure time. 
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If the threshold signal-to-noise ratio is (S/N\h, the threshold 
magnification and the corresponding normalized (by 9e) im- 
pact parameter of the source trajectory for event detection are 
defined respectively as 



and 



«0,th = 



(1.0 



1/2 



FIG. 1 . — Variation of caustic pattern depending on binary separation and example light curves of binary lens events. In each of the left panels, the green dots 
represent the positions of the individual binary lens components, where the heavier lens is toward right. The coordinates are centered at the midpoint between 
the binary lens components. The figure drawn by red curve represents the caustic. The cross mark denotes the effective single lens position, which roughly 
corresponds to the center of each set of caustics. The greyscale represents the signal-to-noise ratio of lensing-induced source flux variation as a function of 
source position. A contour (white curve) is drawn at the level of S/N = 3.0. Presented in the right panels are the light curves of binary lens events resulting 
from the source trajectories marked by arrow in the corresponding left panels. For the case where s = 2.0, we present two light curves and the source trajectories 
responsible for the individual light curves are marked by the corresponding colors. The dark-tone dots on each light curve represent the data points satisfying the 
condition that the source is within the caustic and the signal-to-noise ratio is higher than a threshold value of (S//V)th = 3.0. 

the center of the individual caustics. 10 In Figure Q] the ef- 
fective single lens positions are denoted by cross marks for 
various cases of binary lenses. Once events are produced, we 
estimate the detection efficiency of binary lens events as the 
ratio of the number of events that can be firmly identified to 
be caused by binary lenses out of all produced events. 

Light curves of events are produced under the following ob- 
servational conditions. For the photometry, we assume that 
the instrument can detect 1 photon per second for an / = 24.2 
star following the specification of the Liverpool Telescope. 
The observation is assumed to be carried out such that small- 
exposure images are combined to make a 30 min exposure 
image to obtain a high signal-to-noise ratio while preventing 
saturation in the central bulge region. We assume that 5 such 
combined images are obtained per 24-hour period from the 
observations by using network teles copes. The photometr y 
is done based on difference imaging (lAlard & Luptonl fT998). 
The flux variation is assumed to be measured at an aperture 
that maximizes the signal-to-noise ratio of the measured flux 
variation. Under the assumption of a Gaussian point spread 
function, the optimal apertur e that maximize t he signal-to- 
noise ratio is 6* ap = 0.6736* see dKerins et alj|2006l) . We assume 
that the average seeing is 9 see = 1".0. With the adoption of this 
aperture, the fraction of the source flux within the aperture 
is F(8 < ap )/F tol = 0.715, where F tot is the flux measured at 



(7) 



(8) 



We, therefore, produce events with impact parameters of the 
source trajectory with respect to the effective single lens po- 
sition, Zh,eS, less than the threshold value defined in equa- 
tion ([8]) with an adopted threshold signal-to-noise ratio of 
(S/N)± = 3.0, i.e. 3cr detection. The effective single lens posi- 
tion represents the location of a single lens at which the result- 
ing single lensing light curve best describes the light curve of 
a binary lens event. For example, the light curve of an event 
caused by a close binary lens is well described by the light 
curve of a single lens event caused by a mass equal to the total 
mass of the binary located at the center of mass of the binary. 
In this case, the effective single lens position is the center of 
mass of the binary. For the case of a wide-separation binary, 
on the other hand, the individual lens components behave as if 
they are two independent single lenses located roughly at the 
centers of the individual sets of caustics, i.e. ZL.eff = Zl + Azl- 
In this case, there are two effective lens positions located at 



10 We distinguish wide- separation binaries from close binaries if the caus- 
tic is divided into two pieces. 
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FIG. 2. — Detection efficiency of binary events toward the M31 bulge field 
as a function of the binary separation and mass ratio. The efficiency is for a 
typical pixel-lensing event expected to be detected by a high-cadence pixel- 
lensing survey such as Angstrom. For details about the representative event 
and observational conditions, see § 3. 

# ap = OO. 

We assume the following criteria for the identification of 
binary lens events. The light curves of binary events exhibit 
diverse patterns. In many cases, it is difficult to firmly dis- 
tinguish light curves from those caused by single lenses or 
other types of source variability. Fortunately, a fraction of bi- 
nary lens events involve caustic crossings and the characteris- 
tic features (e.g. spikes during caustic crossings and U-shape 
curve between caustic crossings) in the resulting light curves 
can be used to securely identify not only the lensing-induced 
variability but also the binary nature of lenses. We, there- 
fore, restrict detectable binary events only to caustic-crossing 
events. For the confirmation of caustic crossings, we require 
events should have at least 5 data points with S/N > 3.0 on 
the part of the light curve between caustic crossings. Consid- 
ering that the assumed observation frequency is 5 times per 
day, this requirement is equivalent that source star should stay 
in the caustic at least for a day. We note that the data points 
with S/N > 3.0 do not have to be consecutive. 

In Figure |2j we present the detection efficiency of binary 
pixel-lensing events toward the M31 bulge field as a func- 
tion of the binary separation and mass ratio. From the fig- 
ure, one finds that binaries with separations 0.4 < s < 3.0 can 
be efficiently detected. Especially for binaries with separa- 
tions very close to the Einstein radius, the efficiency is nearly 
100%. One also finds that the dependency of the efficiency on 
the mass ratio is weak and thus detections can be extended to 
low-mass companions. 

4. BINARY-LENS EVENT RATE 

Pixel-lensing events have diverse physical parameters. For 
the estimation of the binary event rate, then, convolution of 
the efficiency estimated for a specific event with distributions 
of the physical lens parameters is required. 

The physical parameters involved wit h primary lens com- 
ponents are based on the simulation of iKerins et alj d2006l) 
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FIG. 3. — Models of binary mass-ratio distribution. 

described in § 3. Once these parameters are set, we t hen se t 
the parameters related to companions. Following lAbtl (1983), 
we model the binary separation as uniformly distributed in 
logs, i.e. /(logs) = const. For the mass ratio distribution, we 
test three models. The first model is based on the assumption 
that the two masses of the binary components are drawn inde- 
pendently from the same mass function as that of single stars. 
This model is the natural result of binary formation process 
where binaries are formed th rough interactions between pro- 
tostellar disks dPringldri989l) or some other form of capture. 
We refer to this model as the 'capture' model. Other possible 
mechanisms of binary formation are fission of a single star 
and fragmentation of a collapsing object. Numerical calcula- 
tions suggest that the former process results in a mass ratio 
distribution peaking at around q = 0.2 dLucvll 19771). while th e 
latter results in more equal masses (Nor man & Wilsonlll978l) . 
We refer to these models as the 'fission' and 'fragmentation' 
models, respectively. We model the mass ratio distribution in 
the fission model as 

f(V) _/54 for q< 0.2, 

\ — 1.25^+1.25 for^>0.2. w 

The distribution of the fragmentation model is modeled as 

f(q) = q. (10) 

In FigurefJ) we present the mass ratio distributions of the three 
tested models. In the capture and fragmentation models, the 
mass of the primary is drawn from the single lens mass func- 
tion and the mass of the companion is determined based on 
the mass ratio derived from the mass ratio distribution. In 
the fission model, on the other hand, we set the total mass of 
the binary to be the same as that of a single lens because a 
single mass is split into two components of a binary in this 
model. Then, the total mass of the binary is drawn from the 
mass function of single lenses and the masses of the individ- 
ual binary components are deduced from the mass ratio. As 
a result, the average mass of the binary lenses in the fission 
model is smaller than those in the capture and fragmentation 
models. 
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FIG. 4. — Distributions of mass ratios (upper panel), normalized separations 
(middle panel), and physical separations (bottom panel) of binary lenses ex- 
pected to be detected by a high-cadence pixel-lensing survey towards M3 1 . 
The three histograms in each panel represent the distributions under different 
models of binary mass-ratio distributions. 



TABLE 1 
Rate of Binary Lens Events 



mass ratio model 



binary event rate 



capture model 
fission model 
fragmentation model 



12.0 /i, events per season 
7.4 /i, events per season 
15.2 /b events per season 



NOTE. — Rates of binary lens events expected 
to be detected by a high-cadence pixel-lensing sur- 
vey toward M31 field under three different mod- 
els of binary mass-ratio distributions. The quantity 
/b represents the fraction of binaries with projected 
separations of 1(T 3 AU < d < 10 3 AU out of all 
lenses. 

In Table [T] we present the expected rates of binary lens 
events under three different models of binary mass ratio dis- 
tribution. We note that the event rate is normalized so that it 
becomes 50 events per season if all lenses are composed of 
a sing le component following the estimation of iKerins et alj 
(2006). We find that the rate of binary lens events is 1\, ~ 7/b 
- 15/b events per season, where ft, is the fraction of binary 
lenses with projected separations of 10~ 3 AU < d < 10 3 AU 
out of all lenses. By adopting a value of ft, = 0.5 and con- 
sidering the ~ 7 -month duration of M3 1 season (from August 
through to February), this rate roughly corresponds to ~ 0.7 
events per month. We find that the dependence of the detec- 
tion rate on the binary mass ratio distribution is such that the 
rate increases as more companions are populated in higher 
mass-ratio region. As a result, the rate becomes bigger in 
the order of the fission, capture, and fragmentation models. 



In Figure [4] we also present the distributions of mass ratios 
and separations (both normalized and physical separations) 
between the components of binary lenses. From the figure, 
we find that binaries to be detected by high-cadence pixel- 
lensing surveys will have mass ratios distributed over a wide 
range of q > 0.1. By contrast, separations of binaries will be 
populated within a narrow range of 1 AU < d < 5 AU. 

5. DISCUSSION AND CONCLUSION 

We estimated the detection rate of binary lens events ex- 
pected from high-cadence pixel-lensing surveys toward M3 1 
such as the Angstrom Project based on detailed simulation 
of events and application of realistic observational condi- 
tions. Under the conservative detection criteria that only 
high signal-to-noise caustic-crossing events with long enough 
durations between caustic crossings can be firmly identified 
as binary lens events, we estimated that the rate would be 
Tb ~ (7- 15)/b(iV/50) per season, where N is the rate of 
stellar pixel-lensing events. We found that detected bina- 
ries would have mass ratios distributed over a wide range of 
q > 0.1 but with separations populated within a narrow range 
of 1 AU < c/ < 5 AU. 

The Angstrom survey is currently commissioning operation 
with a real-time data processing pipeli ne and a web-base d 
transient and microlensing alert system (Da rnlev et al.ll2007l) . 
Realization of the alert system and subsequent follow-up ob- 
servations would greatly increase the observation frequency 
and thus overall event rate. Considering that a significant 
fraction of caustic crossing events are not detected from de- 
tections due to the short duration between caustic crossings, 
increasing the monitoring frequency from follow-up observa- 
tions would be able to dramatically increase the binary event 
rate. 

Follow-up observation of binary lens events would be im- 
portant not only for the increase of the event rate but also 
for the characterization of lens matter. The majority of M3 1 
pixel-lensing events are associated with giant stars with large 
angular radii. The large source size allows resolution of caus- 
tic crossings because the duration of the caustic crossing is 
proportional to the source radius, i.e. 



Af c , 



P* 



(11) 



where = 9±/9e is the normalized source radius and a is 
the incidence angle of the source trajectory with respect to 
the caustic curve. For a typical M31 pixel-lensing event with 
Pi, ~ 0.03 and ?e ~ 10 days, the caustic-crossing duration is 
A? cc ^ 14 hours. Then, if events can be followed up with a 
frequency higher than once every hour, the caustic crossing 
can be resolved. Once the value of p+ is known from the anal- 
ysis of the light curve, one can estimate the Einstein radius 
with additional information of the source radius by 9e = 9*/ p*. 
Since the Einstein radius is related to the physical parameters 
of the lens by the relation in equation (f2j), one can better con- 
strain the nature of the lens. 



CH and B-GP are supported by the grant (C00072) of the 
Korea Research Foundation. Y-BJ and C-UL acknowledge 
the support from Korea Astronomy and Space Science Insti- 
tute. MI acknowledges the support from the grant R01-2006- 
00- 106 10-0 provided by the Basic Science Research Program 
of the Korea Science and Engineering Foundation. EK and 
JPD are supported, respectively, by the Advanced Fellowship 



6 



PIXEL-LENSING DETECTIONS OF M31 BINARIES 



and PhD studentship from the UK Particle Physics and As- 
tronomy Research Council. MAI and RGK acknowledge the 
permanent technical and financial support of the Maidanak 



Observatory operations from the consortium of Korea Uni- 
versities operating under a MOU between the consortium and 
UBAI. 



REFERENCES 



Abt, H. A. 1983, ARA&A, 21, 343 
Afonso, C, et al. 2003, A&A, 400, 951 
Alard, C, & Lupton, R. 1998, ApJ, 503, 325 
Alcock, C, et al. 2000, ApJ, 542, 281 
Alcock, C, et al. 2000, ApJ, 542, 281 
Alcock, C, et al. 2001, ApJ, 557, 1035 
Ansari, R., et al. 1999, A&A, 344, L49 
Baltz, E. A., & Gondolo, P. 2001, ApJ, 559, 41 
Bond, I. A., et al. 2001, MNRAS, 327, 868 
Calchi Novati, S., et al. 2005, A&A, 443, 91 1 
Darnley, M. J., et al. 2007, ApJ, submitted 
de Jong, J. T. A., et al. 2004, A&A, 417, 461 
Di Stefano, R., & Mao, S. 1996, ApJ, 457, 93 
Hamadache, C, et al. 2006, A&A, 454, 185 

Kerins, E., Darnley, M. J., Duke, J. P., Gould, A., Han, C, Jeon, Y.-B.. 
Newsam, A., & Park, B.-G. 2006, MNRAS, 365, 1099 



Lucy, L. B. 1977, AJ, 82, 1013 
Mao, S., & Paczyriski, B. 1991, ApJ, 374, L37 
Norman, M. L., & Wilson, J. R. 1978, ApJ, 224, 497 
Pringle, J. E. 1989, MNRAS, 239, 361 

Riffeser, A., Fliri, J., Bender, R., Seitz, S., & Gossl, C. A. 2003, ApJ, 599, 
L17 

Sumi, T., et al. 2006, ApJ, 636 240 

Tisserand, P., et al. 2006, A&A, submitted I astro-ph/0607207 I 

Uglesich, R. R., Crotts, A. P. S., Baltz, E. A., de Jong J., Boyle, R. P., & 

Corbally, G. J. 2004, ApJ, 612, 877 
Witt, H. J. 1990, A&A, 236, 311 



